Predictive Interpolation of a Video Signal

ABSTRACT

A method of enhancing picture quality of a video signal is disclosed. The method comprises the steps of generating an interpolated base frame image; receiving a first previously decoded difference picture; receiving a second previously decoded difference picture; generating a combined motion compensated difference surface; and generating a temporally interpolated enhanced picture based upon the interpolated base frame image and the combined motion compensated difference surface. A circuit for enhancing picture quality of a video signal is also disclosed. The circuit comprises a base decoder generating a base image of a standard definition picture; a temporal predictive interpolator coupled to the base decoder and generating an interpolated block; and a summing circuit coupled to the temporal predictive interpolator. The summing circuit preferably adds the interpolated block and a difference block.

CLAIM FOR PRIORITY

This application is a divisional application of prior pending U.S.patent application Ser. No. 11/782,915, entitled “PredictiveInterpolation of A Video Signal”, filed Jul. 25, 2007, which is acontinuation of prior U.S. patent application Ser. No. 10/446,347, filedon May 28, 2003, now U.S. Pat. No. 7,386,049, which claims priority toU.S. Provisional Application 60/384,047, entitled “Video InterpolationCoding”, filed on May 29, 2002. All of which are hereby incorporated byreference in their entireties.

RELATED APPLICATIONS

This application relates to U.S. application Ser. No. 10/447,213,entitled “Video Interpolation Coding”, U.S. application Ser. No.10/447,296, entitled “Classifying Image Areas of A Video Signal”, andU.S. application Ser. No. 10/447,216, entitled “Maintaining A pluralityof Codebooks Related to A Video Signal”, each of which beingincorporated herein by reference.

BACKGROUND

Pixonic High Definition (PHD) significantly improves perceptual detailof interpolated digital video signals with the aid of a small amount ofenhancement side information. In its primary application, PHD rendersthe appearance of High Definition Television (HDTV) picture quality froma standard Definition Television (SDTV) coded DVD movie which has beenoptimized, for example, for a variable bit rate average around 6 mbps(megabits-per-second), while the multiplexed enhancement stream averagesapproximately 2 mbps.

In 1953, the NTSC broadcast system added a scalable andbackwards-compatible color sub-carrier signal to then widely deployed525-line black-and-white modulation standard. Newer television receiversthat implemented NTSC were equipped to decode the color enhancementsignal, and then combine it with the older black-and-white componentsignal in order to create a full color signal for display. At the sametime, neither the installed base of older black-and-white televisions,nor the newer black-and-white only televisions designed withforeknowledge of NTSC would need color decoding circuitry, nor would benoticeably affected by the presence of the color sub-carrier in themodulated signal. Other backwards-compatible schemes followed NTSC.

Thirty years later, PAL-Plus (ITU-R BT.1197) added a sub-carrier to theexisting PAL format that carries additional vertical definition forletterboxed video signals. Only a few scalable analog video schemes havebeen deployed, but scalability has been more widely adopted in audiobroadcasting. Like FM radio, the North American MTS stereo (BTSC) audiostandards for television added a sub-carrier to modulate the stereodifference signal, which when matrix converted back to discrete L+Rchannels, could be combined in advanced receivers with the mono carrierto provide stereo audio.

In most cases, greater spectral efficiency would have resulted if theencoding and modulation schemes had been replaced with state-of-the-artmethods of the time that provided the same features as the scalableschemes. However, each new incompatible approach would have displacedthe installed base of receiving equipment, or required spectruminefficient simulcasting. Only radical changes in technology, such asthe transition from analog to digital broadcast television, haveprompted simultaneous broadcasting (“simulcasting”) of related content,or outright replacement of older equipment.

Prior attempts to divide a compressed video signal into concurrentscalable signals containing a base and at least one enhancement layerhave been under development since the 1980's. However, unlike analog, nodigital scalable scheme has been deployed in commercial practice,largely due to the difficulties and overheads created by the scalabledigital signals. The key reason perhaps is found is in the very naturein which the respective analog and digital consumer distribution signalsare encoded: analog spectra have regular periods of activity (orinactivity) where the signal can be cleanly partitioned, while digitalcompressed signals have high entropy and irregular time periods thatcontent is modulated.

Analog signals contain high degree of redundancy, owing to theirintended memory-less receiver design, and can therefore be efficientlysliced into concurrent streams along arbitrary boundaries within thesignal structure. Consumer digital video distribution streams such asDVD, ATSC, DVB, Open Cable, etc., however apply the full toolset ofMPEG-2 for the coded video representation, removing most of theaccessible redundancy within the signal, thereby creating highlyvariable, long-term coding dependencies within the coded signal. Thisleaves fewer cleaner dividing points for scalability.

The sequence structure of different MPEG picture coding types (I, P, B)has a built-in form of temporal scalability, in that the B pictures canbe dropped with no consequence to other pictures in the sequence. Thisis possible due to the rule that no other pictures are dependently codedupon any B picture. However, the instantaneous coded bitrate of picturesvaries significantly from one picture to another, so temporal scalablebenefits of discrete streams is not provided by a single MPEG bitstreamwith B-pictures.

The size of each coded picture is usually related to the content, orrate of change of content in the case of temporally predicted areas ofthe picture. Scalable streams modulated on discrete carriers, for thepurposes of improved broadcast transmission robustness, aretraditionally designed for constant payload rates, especially when asingle large video signal, such as HDTV, occupies the channel. VariableBit Rate (VBR) streams provide in practice 20% more efficient bitutilization that especially benefits a statistical multiplex ofbitstreams.

Although digital coded video for consumer distribution is only a recentdevelopment, and the distribution mediums are undergoing rapidevolution, such as higher density disks, improved modems, etc., scalableschemes may bridge the transition period between formats.

The Digital Versatile Disc (DVD), a.k.a. “Digital Video Disc,” format isdivided into separate physical, file systems, and presentation contentspecifications. The physical and file formats (Micro-UDF) are common toall applications of DVD (video, audio only, computer file). Video andaudio-only have their respective payload specifications that define thedifferent data types that consume the DVD storage volume.

The video application applies MPEG-2 Packetized Elementary Streams (PES)to multiplex at least three compulsory data types. The compulsory streamtypes required by DVD Video are: MPEG-2 Main Profile@Main Level(standard definition only) for the compressed video representation;Dolby AC-3 for compressed audio; a graphic overlay (sub-picture) format;and navigation information to support random access and other trick playmodes. Optional audio formats include: raw PCM; DTS; and MPEG-1 LayerII. Because elementary streams are encapsulated in packets, and asystems demultiplexer with buffering is well defined, it is possible forarbitrary streams types to be added in the future, without adverselyaffecting older players. It is the role of the systems demultiplexer topass only relevant packets to each data type specific decoder.

Future supplementary stream types envisioned include “3D” stereo vision,metadata for advanced navigation, additional surround-sound ormultilingual audio channels, interactive data, and additional videostreams (for supporting alternate camera angles) that employ moreefficient, newer generation video compression tools.

Two major means exist for multiplexing supplementary data, such asenhancement stream information of this invention, in abackwards-compatible manner. These means are not only common to DVD, butmany other storage mediums and transmission types including D-VHS,Direct Broadcast Satellite (DBS), digital terrestrial television (ATSC &DVB-T), Open Cable, among others. As the first common means, the systemsstream layer multiplex described above is the most robust solution sincethe systems demultiplexer, which comprises a parser and buffer, iscapable of processing streams at highly variable rates withoutconsequence to other stream types multiplexed within the same systemsstream. Further, the header of these system packets carry a uniqueRegistered ID (RID) that, provided they are properly observed by thecommon users of the systems language, uniquely identify the stream typeso that no other data type could be confused for another, includingthose types defined in future. SMPTE-RA is such an organization chargedwith the responsibility of tracking the RID values.

The other, second means to transport supplementary data, such asenhancement data of the invention, is to embed such data within theelementary video stream. The specific such mechanisms available toMPEG-1 and MPEG-2 include user_data( ) extension start codes, reservedstart codes. Other coding languages also have their own means ofembedding such information within the video bitstream. These mechanismshave been traditionally employed to carry low-bandwidth data such asclosed captioning and teletext. Embedded extensions provides a simple,automatic means of associating the supplementary data with the intendedpicture the supplementary data relates to since these embedded transportmechanisms exist within the data structure of the correspondingcompressed video frame. Thus, if a segment of enhancement data is foundwithin a particular coded picture, then it is straight-forward for asemantic rule to assume that such data relates to the coded picture withwhich the data was embedded. Also, there is no recognized registrationauthority for these embedded extensions, and thus collisions betweenusers of such mechanisms can arise, and second that the supplementarydata must be kept to a minimum data rate. ATSC and DVD have madeattempts to create unique bit patterns that essentially serve as theheaders and identifiers of these extensions, and register the ID's, butit is not always possible to take a DVD bitstream and have it translatedirectly to an ATSC stream.

Any future data stream or stream type therefore should have a uniquestream identifier registered with, for example, SMPTE-RA, ATSC, DVD,DVB, OpenCable, etc. The DVD author may then create a PacketizedElementary Stream with one or more elementary streams of the this type.

Although the sample dimensions of the standard definition format definedby the DVD video specification are limited to 720×480 and 720×576 (NTSCand PAL formats, respectively), the actual content of samples may besignificantly less due to a variety of reasons.

The foremost reason is the “Kell Factor,” which effectively limits thevertical content to approximately somewhere between ⅔ and % response.Interlaced displays have a perceived vertical rendering limit between300 and 400 vertical lines out of a total possible 480 lines of content.DVD video titles are targeted primarily towards traditional 480i or 576idisplays associated with respective NTSC and PAL receivers, rather thanmore recent 480p or computer monitors that are inherently progressive(the meaning of “p” in 480p). A detailed description of the Kell Factorcan be found in the books “Television Engineering Handbook” by Wilkonsonet al, and “Color Spaces” by Charles Poynton. A vertical reduction ofcontent is also a certain measure to avoid the interlace flicker problemimplied by the Kell Factor. Several stages, such as “film-to-tape”transfer, can reduce content detail. Interlace cameras often employlenses with an intentional vertical low-pass filter.

Other, economical reasons favor moderate content reduction.Pre-processing stages, especially low-pass filtering, prior to the MPEGvideo encoder can reduce the amount of detail that would need to beprescribed by the video bitstream. Assuming, the vertical content isalready filtered for anti-flicker (Kell Factor), filtering along thehorizontal direction can further lower the average rate of the codedbitstream by a factor approximately proportional to the strength of thefiltering. A 135 minute long movie would have an average bitrate of 4mbps if it were to consume the full payload of a single-sided,single-layer DVD (volume of 4.7 billion bytes). However, encoding of720×480 interlace signals have been shown to require sustained bitratesas high as 7 or 8 mbps to achieve transparent orjust-noticeable-difference (JND) quality, even with a well-designedencoder. Without pre-filtering, a 4 mbps DVD movie would likelyotherwise exhibit significant visible compression artifacts. Themeasured spectral content of many DVD tiles is effectively less than 500horizontal lines wide (out of 720), and thus the total product (assuming350 vertical lines) is only approximately half of the potentialinformation that can be expressed in a 720×480 sample lattice. It is notsurprising then that such content can fit into half the bitrate impliedat least superficially by the sample lattice dimensions.

The impact of this softening is minimized by the fact that most 480itelevision monitors are not capable of rendering details within theNyquist limits of 720×480. The displays are likely optimized for aneffective resolution of 500×350 or worse. Potentially, anti-flickerfilters, as commonly found in computer-to-television format converters,could be included in every DVD decoder or player box, thus allowing true480 “p” content to be encoded on all DVD video discs. Such a usefulfeature was neither given as a mandate nor suggested as an option in theoriginal DVD video specification. The DVD format was essentially seen asa means to deliver the best standard definition signals of the time toconsumers.

Prior art interpolation methods can interpolate a standard definitionvideo signal to, for example, a high definition display, but do not addor restore content beyond the limitations of the standard-definitionsampling lattice. Prior art methods include, from simplest to mostcomplex: sample replication (“zero order hold”), bi-linearinterpolation, poly-phase filters, spline fitting, POCS (Projection onConvex Sets), and Bayesian estimation. Inter-frame methods such assuper-resolution attempt to fuse sub-pixel (or “sub-sample”) detail thathas been scattered over several pictures by aliasing and other diffusionmethods, and can in fact restore definition above the Nyquist limitimplied by the standard definition sampling lattice. However suchschemes are computationally expensive, non-linear, and do not alwaysyield consistent quality gains frame-to-frame.

The essential advantage of a high-resolution representation is that itis able to convey more of the actual detail of a given content than alow-resolution representation. The motivation of proving more detail tothe viewer is that it improves enjoyment of the content, such as thequality difference experienced by viewers between the VHS and DVDformats.

High Definition Television (HDTV) signal encoding formats are a directattempt to bring truly improved definition, and detail, inexpensively toconsumers. Modern HDTV formats range from 480p up to 1080p. This rangeimplies that content rendered at such resolutions has anywhere from twoto six times the definition as the traditional, and usually diluted,standard definition content. The encoded bitrate would also becorrespondingly two to six times higher. Such an increased bitrate wouldnot fit onto modem DVD volumes with the modern MPEG-2 video codinglanguage. Modern DVDs already utilize both layers, and have only enoughroom left over for a few short extras such as documentaries and movietrailers.

Either the compression method or the storage capacity of the disc wouldhave to improve to match as the increase in definition and correspondingbitrate of HDTV. Fortunately both storage and coding gains have beenrealized. For example, H.264 (a.k.a. MPEG-4 Part 10 “Advanced VideoCoder”) has provided a nominal 2× gain in coding efficiency over MPEG-2.Meanwhile, blue-laser recording has increased disc storage capacity byat least 3× over the original red-laser DVD physical format. The minimalcombined coding and physical storage gain factor of 6:1 means that it ispossible to place an entire HDTV movie on a single-sided, single-layerdisc, with room to spare.

A high-definition format signal can be expressed independently(simulcast) or dependently (scalable) with respect to astandard-definition signal. The simulcast method codes the standarddefinition and high definition versions of the content as if they wereseparate, unrelated streams. Streams that are entirely independent ofeach other may be multiplexed together, or transmitted or stored onseparate mediums, carriers, and other means of delivery. The scalableapproach requires the base stream (standard definition) to be firstdecoded, usually one frame at a time, by the receiver, and then theenhancement stream (which generally contains the difference informationbetween the high definition and standard definition signals) to bedecoded and combined with the frame. This may be done piecewise, as forexample, each area of the base picture may be decoded just in time priorto the addition of the enhancement data. Many implementation schedulesbetween the base and enhancement steps are possible.

The simulcast approach is cleaner, and can be more efficient thanenhancement coding if the tools and bitrate ratios between the two arenot tuned properly. Empirical data suggests that some balance of ratesshould exist between the base and enhancement layers in order to achieveoptimized utilization of bits. Thus if a data rate is required toachieve some picture quality for the base layer established by theinstalled base of DVD players, for example, then the enhancement layermay require significant more bits in order to achieve a substantialimprovement in definition.

In order to lower the bitrate of the enhancement layer, several trickscan be applied that would not noticeably impact quality. For example,the frequency of intra pictures can be decreased, but at the tradeoff ofreduced robustness to errors, greater IDCT drift accumulation, andreduced random access frequency.

Previous scalable coding solutions have not been deployed in main-streamconsumer delivery mediums, although some forms of scalability have beensuccessfully applied to Internet streaming. With the exception oftemporal scalability (FIG. 2 e) that is inherently built-in all MPEGbitstreams that utilize B-frames, the spatial scalable scheme (FIG. 2d), SNR scalable (FIG. 2 c) and Data Partitioning schemes documented inthe MPEG-2 standard have all incurred a coding efficiency penaltyrendering scalable coding efficiency little better, or even worse, thanthe total bandwidth consumed by the simulcast approach (FIG. 2 b). Thereasons behind the penalties have not been adequately documented, butsome of the known factors include: excessive block syntax overheadincurred when describing small-enhancements, and re-circulation ofquantization noise between the base and enhancement layers.

FIG. 2 a establishes the basic template where, in subsequent figures,the different scalable coding approaches most fundamentally differ intheir structure and partitioning. Bitstream Processing (BP) 2010includes those traditional serially dependent operations that have avarying density of data and hence variable complexity per coding unit,such as stream parsing, Variable Length Decoding (VLD), Run-LengthDecoding (RLD), header decoding. Inverse Quantization (IQ) is sometimesplaced in the BP category if only the non-zero transform coefficientsare processed rather applying a matrix operation upon all coefficients.Digital signal processing (DSP) 2020 operations however tend to beparallelizable (e.g. SIMD scalable), and have regular operations andcomplexity. DSP includes IDCT (Inverse Discrete Cosine Transform) andMCP (Motion Compensated Prediction). Reconstructed blocks 2025 arestored 2030 for later display processing (4:2:0 to 4:2:2 conversion,image scaling, field and frame repeats) 2040, and to serve as referencefor prediction 2031. From the bitstream 2005, the BP 2010 producesIntermediate decoded bitstream 2015 comprising arrays of transformcoefficients, reconstructed motion vectors, and other directives thatwhen combined and processed through DSP produce the reconstructed signal2025.

FIG. 2 b demonstrates the “simulcast” case of two independent streamsand decoders that optionally, through multiplexer 2136, feed the seconddisplay processor 2140. The most typical application fitting the FIG. 2b paradigm is a first decoder system for SDTV, and a second decodersystem for HDTV. Notably, the second decoder's BP 2110 and DSP 2120stages do not depend upon state from the first decoder.

The scalable schemes are best distinguished by what processing stagesand intermediate data they relate with the base layer. The relationpoint is primarily application-driven. FIG. 2 c illustrates frequencylayering, where the relation point occurs at the symbol stages prior toDSP. (symbols are an alternate name for bitstream elements). In blockbased transform coding paradigms, the symbol stream is predominately inthe frequency domain, hence frequency layering. The enhancedintermediate decoded symbols 2215 combined with the intermediate decodedbase symbols 2015 creates a third intermediate symbol stream 2217 thatis forward-compatible decodable, in this example, by the base layer DSPdecoder 2220. The combined stream appears as an ordinary base layerstream with increased properties (bitrate, frame rate, etc.) over thebase stream 2005. Alternatively, the enhanced DSP decoder could havetools not present in the base layer decoder DSP, and 2217 depending onthe tools combination and performance level, therefore only bebackward-compatible (assuming the enhanced DSP is a superset of the baseDSP). SNR scalability and Data partitioning are two known cases offrequency layering that produce forward-compatible intermediate datastreams 2217 decodable by base layer DSP stages 2020. Frequency layeringis generally chosen for robustness over communications mediums.

In a forward-compatible application example of frequency layering,detailed frequency coefficients that could be added directly to the DCTcoefficient block would be encoded in the enhancement stream, and added2216 to the coefficients 2015 to produce a higher fidelity reconstructedsignal 2225. The combined stream 2217 resembles a plausible base layerbitstream coded at a higher rate, hence the forward compatibledesignation. Alternatively, a backward-compatible example would be anenhancement stream that inserted extra chrominance blocks into thebitstream in a format only decodable by the enhanced DSP decoder. Theoriginal Progressive JPEG mode and the more recent JPEG2000 are examplesof frequency layering.

Spatial scalability falls into the second major scalable codingcategory, spatial layering, whose basic decoding architecture is shownin FIG. 2 d. The spatial scalability paradigm exploits the base layerspatial-domain reconstruction 2025 as a predictor for the enhancedreconstruction signal 2327, much like previously reconstructed picturesserve as reference 2031 for future pictures (only the reference picturesin this example are, as an intermediate step, scaled in resolution). Atypical application would have the base layer contain a standarddefinition (SDTV) signal, while the enhancement layer would encode thedifference between the scaled high definition (HDTV) and standarddefinition reconstruction 2025 scaled to match the lattice of 2325.

Spatial layering is generally chosen for scaled decoder complexity, butalso serves to improve robustness over communications mediums when thesmaller base layer bitstream is better protected against errors in thecommunications channel or storage medium.

A third scalability category is temporal layering, where the base layerproduces a discrete set of frames, and an enhancement layer addsadditional frames that can be multiplexed in between the base layerframes. An example application is a base layer bitstream consisting ofonly I and P pictures could be decoded independently of an enhancementstream containing only B-pictures, while the B-pictures would bedependent upon the base layer reconstruction, as the I and P framereconstructions would serve as forward and backward MCP (MotionCompensated Prediction) references. Another application is stereovision, where the base layer provides the left eye frames, and theenhancement layer predicts the right eye frames from the left eyeframes, with additional correction (enhancement) to code the left-rightdifference.

Enhancement methods that do not employ side information or anysignificant enhancement layer stream are applied by default in theconversion of SDTV to HDTV. Interpolation, through scaling andsharpening, a standard definition (SDTV) signal to a high definition(HDTV) signal is a method to simulate high definition content, necessaryto display SDTV on a high definition monitor. Although the result willnot look as good as genuine HDTV content, certain scaling orinterpolation algorithms do a much better job than others, as somealgorithms better model the differences between a HDTV and SDTVrepresentation of the same content. Edges and textures can be carefullysharpened to provide some of the appearance of HDTV, but will at thesame time look artificial since the interpolation algorithm will notsufficiently estimate the true HDTV from the content. Plausible detailpatterns can be substituted, but may also retain a synthetic look uponclose examination.

Many methods falling under the genre of superresolution can partiallyrestore HDTV detail from an SDTV signal under special circumstances,although to do so requires careful and complex motion compensatedinterpolation since the gain is realized by solving for detail that havebeen mixed over several pictures through iterative mathematicaloperations. Superresolution tools require sub-pixel motion compensatedprecision, similar to that found in newer video coders, and withprocessing at sub-pixel granularity rather than whole blocks. Thus,instead of one motion vector for every 8×8 block (every 64 pixels),there would be one to four motion vectors generated by thesuperresolution restoration algorithm at the receiver for everyhigh-definition pixel. Optimization techniques can reduce thiscomplexity, but the end complexity would nonetheless exceed the combineddecoding and post-processing complexity of the most advanced consumervideo systems. In an effort to improve stability of the restored image,and reduce implementation costs, several approaches have beeninvestigated by researchers to restore high resolution from acombination of a lower resolution image and side information or explicitknowledge available only to the encoder.

Gersho's 1990 publication “non-linear VQ interpolation.” [Gersho90]first proposes to interpolate lower resolution still images by means ofVector Quantization (VQ) codebooks (2410 and 2516) trained on theiroriginal higher resolution image counterparts. Prior interpolationmethods, such as multi-tap polyphase filter banks, generate theinterpolated image sample-by-sample (or point-wise) where data is fittedto a model of the interpolated signal through convolution with curvesderived from the model. The model is typically a sine function. Gersho'sinterpolation procedure (FIG. 20 closely resembles block coding, wherethe picture (example shown in FIG. 7 e) is divided into a grid of inputblocks similar to the grid 7411. Each block (whose relationship to thegrid 7411 is demonstrated by block 7431) in signal 2506 may be processedindependently of other blocks within the same picture. The mapping stage2504 models some form of distortion such as sub-sampling of the originalsignal 2502 to the input signal 2506. It is the goal of the Gersho90interpolator that the reconstructed block 2518 best approximates theoriginal block 2502 given the information available in the receiver,namely, input block 2506 and previously derived codebooks 2510 and 2516.Input block 2506 is matched to a best-fit entry within a first codebook2510. FIG. 2 g adapts the mapping stage 2604 as a combination ofdecimation followed by the MPEG encode-decode process, the focus of thisdisclosure's application. Specifically, the mapping stage is theconversion of an HDTV signal to an SDTV signal (via sub-sampling ordecimation) that is then MPEG encoded. While the classic VQ picturecoder transmits codebook indices to the receiver, in the nonlinear VQinterpolation application (FIG. 2 f through 2 i), the first index 2512of the matching codebook entry in 2510 serves as the index of acorresponding entry in a second codebook 2516. “Super-resolution” isachieved in that the second codebook contains detail exceeding thedetail of the input blocks 2506. Gersho90 is targeted for theapplication of image restoration, operating in a receiver that is giventhe distorted image and codebooks 2510, 2516, 2610, and 2616 trained oncontent 2502 available only at the transmitter.

Gersho's non-linear VQ interpolation method is applied for imagerestoration, and therefore places the codebook search matching and indexcalculation routine at the receiver. In contrast, the typicalapplications of VQ are for compression systems whose search routine isat the transmitter where indices and the codebooks are generated andtransmitted to the receiver. The receiver then uses the transmittedelements to reconstruct the encoded images. While in the Gersho90design, the index generator 2008 is the receiver, the codebook generatorstill resides at the transmitter, where the higher resolution sourcecontent 2002 upon which C* (2016, 2116) is trained, is available.

The principal step of Non-linear Interpolative Vector Quantization forImage Restoration described by [Sheppard97], over the [Gersho90] paperthat it builds upon, is the substitution of the first VQ stage(2508,2608) with a block waveform coder comprising a Discrete CosineTransform 2904 and transform coefficient Quantization stage 2908. Thequantized coefficients are packed 2912 to form the index 2914 applied tothe second codebook 2716, 2812. Thus, a frequency domain codebook iscreated rather than the traditional, spatial domain VQ codebook. Thesignificance of this step is many-fold. First, the codebook searchroutine is reduced to negligible complexity thanks to the combination ofDCT, quantization, and packing stages (2904, 2908, 2912 respectively)that collectively calculate the second codebook index 2712 directly froma combination of quantized DCT coefficients 2906 within the same block2902. Prior methods, such as Gersho90, generated the index through acomprehensive spatial domain match tests (similar to the process in5400) of many codebook entries (similar to 5140) to find the best match,where the index 2712 of the best match serves as the index sought by thesearch routine.

Sheppard further overlaps each input block by a pre-determined number ofsamples. Thus, a window of samples is formed around the projected areato be interpolated, and the input window steps through the picture at anumber of samples smaller than the dimensions of the input block.Alternatively, in a non-overlapping arrangement, the projected and inputblock dimensions and step increments would be identical. An overlappingarrangement induces a smoothing constraint, resulting in a more accuratemapping of input samples to their output interpolated counterparts. Thisleads to fewer discontinuities and other artifacts in the resultinginterpolated image. However, the greater the overlap, the moreprocessing work must be done in order to scale an image of a given size.For example, in a combination of a 4×4 process block overlapping a 2×2input block, sixteen samples are processed for every four samples thatare interpolated. This is a 4:1 ratio of process bandwidth to inputwork. In a non-overlapping arrangement, sixteen samples (in a 4×4 block)are produced for every sixteen input samples. The overlapping examplegiven here requires four times as much work per average output sample asthe non-overlapping case.

Although the DCT method by Sheppard et al does permit larger codebooksthan the NLIVQ methods of Gersho et al, it does not address the cost anddesign of sending such codebooks to a receiver over a communications orstorage medium. The application is a “closed circuit” system, withvirtually unlimited resources, for restoring images of similarresolution. Thus, an improved system that is designed specifically forentropy-constrained, real-time transmission and can scale across imageresolutions is needed.

DVD

DVD is the first inexpensive medium to deliver to main stream consumersnearly the full quality potential of SDTV. Although a rigid definitionof SDTV quality does not exist, the modem definition has settled on“D-1” video—the first recording format to adopt CCM 601 parameters. SDTVquality has evolved significantly since the first widespreadintroduction of television in the 1940's, spawning many shades ofquality that co-exist today.

In the late 1970's, the first popular consumer distribution format, VHSand Betamax tape, established the most common denominator for standarddefinition with approximately 250 horizontal luminance lines and asignal-to-noise ratio (SNR) in the lower to mid 40's dB range. Earlytelevision broadcasts had similar definition. In the 1980's, televisionmonitors, analog laserdiscs, Super-VHS and the S-Video connector offeredconsumers improved SD video signals with up to 425 horizontal lines andSNR as high as 50 dB, exceeding the 330horizontal-line-per-picture-height limit of the broadcast NTSC signalformat today.

Starting in 1982, professional video engineering organizationscollaborated on the creation of the CCIR 601 discrete signalrepresentation standard for the exchange of digital signals betweenstudio equipment. Although it is only one set of parameters among manypossible choices, CCIR 601 effectively established the upper limit forstandard definition at 540 horizontal lines per picture height (on a 4:3aspect ratio monitor). Applications such as DVD later diluted the samepixel grid to cover a one third wider screen area. Thus the horizontaldensity on 16:9 anamorphic DVD titles is one third less than standard4:3 “pan & scan” titles. The CCIR. 601 rectangular grid sample latticewas defined as 720 samples per line, with approximately 480 lines perframe at the 30 Hz frame rate most associated with NTSC, and 576 linesat the 25 Hz frame rate of PAL and SECAM. Horizontal line density iscalculated as (aspect ratio)*(total lines per picture width). For a 4:3aspect ratio, the yield is therefore ((4/3)*(720))=540 lines.

Although technically a signal format, CCIR 601 cultivated its ownconnotation as the ultimate watermark of “studio quality.” By the late1990's, CCIR 601 parameters were ushered to consumers by the ubiquitousMPEG-2 video standard operating mode, specifically designated “MainProfile@Main Level or “MP@ML”. MPEG-2 MP@ML was adopted as the exclusiveoperating point by products such as DVD, DBS satellite, and digitalcable TV. While the sample dimensions of DVD may be fixed to 720×480(“NTSC”) and 720×576 (“PAL”), the familiar variables such as bitrate(bandwidth), content, and encoder quality very much remain dynamic, andup to the discretion of the content author.

Concurrent to the end of the SDTV evolution, HDTV started from almostits beginning as a handful of digital formats. SMPTE 274M has becomeHDTV's ubiquitous analogy for to SDTV's CCIR 601. With 1920samples-per-line by 1080 lines per frame, and a 16:9 aspect aspectratio—one third wider than the 4:3 ratio of SDTV SMPTE 274M meets thecanonical requirement that HD be capable of rendering twice thehorizontal and vertical detail of SDTV. The second HDTV format, SMPTE296M, has image dimensions of 1280×720 samples.

Until all programming is delivered in an HDTV format, there will be aneed to convert SDTV signals to fit on HDTV displays. SDTV legacycontent may also circulate indefinitely. In order to be displayed on atraditional HDTV display, SDTV signals from sources such as broadcast,VHS, laserdisc, and DVD need to first be up-converted to HDTV. Classicpicture scaling interpolation methods, such as many-tap FIR poly-phasefilters, have been regarded as the state of the art in practicalinterpolation methods. However, the interpolated SD signal will still belimited to the detail prescribed in the original SD signal, regardlessof the sample density or number of lines of the HD display. InterpolatedSD images will often appear blurry compared to their true HDcounterparts, and if the interpolated SD images are sharpened, they maysimulate some aspect of HD at the risk looking too synthetic.

One reason for SD content looking better on HD displays comes from thefact that most display devices are incapable of rendering the fulldetail potential of the signal format they operate upon as input. The HDdisplay has the advantage that details within the SD image that were toofine or subtle to be sufficiently resolved by a SD display can becomemuch more visible when scaled up on the HD display. Early on, however,the interpolation processing and HD display will reach a point ofdiminishing returns with the quality and detail that can be renderedfrom an SD signal. In the end, information must be added to the SDsignal in order to render true detail beyond the native limits of the SDformat. Several enhancement schemes, such as the Spatial Scalable codersof MPEG-2, have been attempted to meet this goal, but none have beendeployed in commercial practice due to serious shortcomings.

Enhancement methods are sensitive to the quality of the base layersignal that they build upon. To optimize the end quality, a balance inbitrate and quality must be struck between the base layer andenhancement layer reconstructions. The enhancement layer should notalways spend bits correcting deficiencies of the base layer, while atthe same time the base layer should not stray too close to its own pointof diminishing returns.

SUMMARY

FIG. 1 a shows the conceptual performance of the invention when used asan enhancement coder in conjunction with an MPEG-2 base layer. Theperceived quality level Q2 achieved with the PHD/MPEG-2 combination atrate R2 is greater than the quality that would be reached using onlyMPEG-2 at the same rate R2. In this figure, MPEG expresses quality up toa natural stopping point, where PHD picks up and carries it further at afaster rate (denoted with a higher Q/R slope). The figure expresses thatthere is a natural dividing point between MPEG-2 and PHD that leads toan overall optimal quality.

While DVD video may be the first popular consumer format to reach thelimits of standard definition, artifacts may still be occasionallyvisible, even on the best coded discs. Those skilled in the art of videocoding are familiar with empirical measures that an MPEG-2 videobitstream can sustain up to 10 million bits per second at transparentquality levels when approximating a CCIR 601 rate standard definitionvideo signal containing complex scenes. Sophisticated pre-processingsteps can be carefully applied to reduce the content of the signal inareas or time periods that will not be very well perceived, andtherefore reduce coded bitrate for those areas, and/or remove datapatterns that would not map to a concise description with the MPEG-2video coding language. Removal of noise, temporal jitter, and film graincan also help reduce bitrate. Human-assisted coding of difficult scenesis used to make decisions on areas or periods that fail encoderanalysis. However, even with these and other optimization steps, theaverage bitrate will, for film content coded at the quality limits ofSDTV, be on the order of 6 to 7 mbps. The reference DVD system, definedby the DVD Forum members and documented in the DVD specification,requires that the DVD player transport and multiplexing mechanism shallindefinitely sustain video rates as high as 9.5 mbps.

Therefore to bridge the transition between the modern DVD standarddefinition format, and any new high definition format that employs acombination of new coding methods and new storage mediums (which are notbackwards compatible with older means), an improved method ofenhancement coding is needed.

The interpolation error signal is the difference between theinterpolated signal and the original signal that the interpolation isattempting to estimate or predict. The interpolation error typically hashigh concentration of energy along edges of objects, since the edges aremost difficult to model accurately with prediction. PHD includes toolsfor the efficient coding of the most perceptible detail within theinterpolation error signal that represents information lost, forexample, during the filtering conversion from the original HD signal tothe base layer signal.

PHD efficiently exploits the base layer video information alreadyavailable to the receiver, thereby minimizing the amount of enhancementinformation to be sent. Two principal tools are employed to this end:the classifier, and the predictive interpolator. In a specific instanceof the preferred embodiment, classification is applied to the base layerto select sub-tables of a codebook that contains a collection ofadditive detail block patterns activated by the coded enhancementstream. The overall algorithm is conceptualized in FIG. 1 b through theillustration of data at various stages of transformation as data passesthrough the PHD decoder.

The preferred instance of the toolset resembles a block-based videocoding language. Difference blocks are first sent within the enhancementbitstream to improve or correct the accuracy of the predicted image.Then, individual blocks are applied to interpolated areas. Small blocksizes, such as the preferred embodiment's 4×4 base layer classificationblock size, offer a reasonable tradeoff between bitrate, implementationcomplexity, and approximation of picture features and contours. Each 4×4area in the base layer image has a corresponding 8×8 area in theinterpolated image.

The PHD decoder analyzes the base layer data, through for example thepreferred classification methods, and adds enhancement data to theinterpolated signal. Many stages of the enhancement process are alsoguided by analysis conducted on the base layer reconstruction. Forexample, flat background areas that are determined unworthy ofenhancement by the base layer analyzer do not incur the overhead ofsignaling in the enhancement stream of how those areas should betreated.

To demonstrate the power of the classification tool, FIG. 1 c shows asmall codebook 1210 of image patterns before and after partitioning byclassification. Codevectors are sorted by their base patterns in theleft column 1210, and then are grouped into the right boxes (1220, 1222,1224, 1226) according to the base pattern common to each cluster ofcodevectors. The simplified example has four codevectors per each of thefour classes. After clustering, the address space 1212 is effectivelycut in half, resulting in a 2-bit index 1221—half the size of theoriginal 4-bit index 1212 (shown along the left column) needed touniquely address each codevector. The first two prefix bits of theoriginal 4-bit index are effectively derived from the base layeranalyzer.

To demonstrate the application of the classifier, FIG. 1 d shows the setof classes for a simple picture with one foreground object (tree) andseveral background areas (sky, mountains, and grass). Each block isassigned a class number in FIG. 1 d, and a separate sub-table codevectorindex in FIG. 1 e. The object outlines in FIG. 1 e illustrate the highpass signal of the solid objects in FIG. 1 d. The high pass, or“difference” signal, is effectively coded with the blocks in thecodebook table.

Any distinct pattern or set of attributes that can be derived from thebase layer, through a combination of operations and analytical stages,and has commonality among a sufficient number of codevectors, can serveas a class. The larger the number of codevectors that share commonattributes (such as the example base patterns in FIG. 1 c), the greaterthe reduction of the global address space of the codebook and hencesmaller the codevector indices that need to be transmitted to the PHDdecoder. In other words, the amount of information that nominally needbe sent can first be reduced by partially deriving whatever informationpossible in the receiver.

Classification also forces unimportant codevectors that do not stronglyfall into any class to merge with like codevectors.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 a is a block diagram showing the performance of the invention.

FIG. 1 b is a block diagram showing the transformation of data as itpasses through a decoder according to the present invention.

FIG. 1 c shows a codebook of image patterns before and afterpartitioning by classification.

FIG. 1 d shows a set of classes for one picture according to the presentinvention.

FIG. 1 e shows a sub-table codevector index according to the presentinvention.

FIG. 2 a shows a block diagram of single non-scalable stream accordingto the present invention.

FIG. 2 b shows a block diagram of two independent streams according tothe present invention.

FIG. 2 c is a block diagram showing frequency layer according to thepresent invention.

FIG. 2 d is a block diagram showing spatial scalability according to thepresent invention.

FIG. 2 e a block diagram showing temporal scalability according to thepresent invention.

FIG. 2 f is a block diagram showing a Gersho interpolation procedure.

FIG. 2 g is a block diagram showing a mapping stage having a combinationof decimation followed by an MPEG encode/decode process according to thepresent invention.

FIG. 2 h is a block diagram showing non-linear interpolation vectorquantization according to the present invention.

FIG. 2 i is a block diagram showing non-linear interpolation vectorquantization of MPEG encoded video.

FIG. 2 j is a block diagram showing index generation steps.

FIGS. 3 a and 3 b are block diagrams showing the fundamental stages of aclassifier according to the present invention.

FIG. 3 c are various blocks and block codings.

FIG. 3 d is a block diagram showing the fundamental stages of aclassifier according to an alternate embodiment of the presentinvention.

FIG. 3 e is a flow chart showing the classification process according tothe present invention.

FIG. 3 f is a block diagram of a state machine according to the presentinvention.

FIG. 4 a is a block diagram showing a conventional spatial scalableenhancement architecture.

FIG. 4 b is a block diagram showing stages of video coding according tothe present invention.

FIG. 4 c is a conventional decoder.

FIG. 4 d is another conventional decoder.

FIG. 4 e is another conventional decoder.

FIG. 4 f is a decoder.

FIG. 4 g is another conventional decoder.

FIG. 5 a is a block diagram of a real-time process stage of anenhancement process according to the present invention.

FIG. 5 b is a block diagram showing databases maintained by an encoderaccording to the present invention.

FIG. 5 c is a block diagram showing look ahead stages of an enhancementencoder according to the present invention.

FIG. 5 d is a block diagram showing a pro-classification stage accordingto the present invention.

FIG. 5 e is a block diagram showing a circuit for authoring figuresaccording to the present invention.

FIG. 5 f.1 is a block diagram showing conventional DVD authorizing.

FIG. 5 f.2 is a block diagram showing DVD authorizing according to thepresent invention.

FIG. 5 g is a block diagram showing storage prior to multiplexing a discrecord.

FIG. 6 a is a block diagram showing stages within the predictionfunction according to the present invention.

FIG. 6 b is a block diagram showing the generation of an enhancedpicture.

FIG. 6 c is a functional block diagram of a circuit for generatingenhanced pictures according to the present invention.

FIG. 6 d is a block diagram of a circuit for generating enhancedpictures according to the present invention.

FIG. 7 a is a strip diagram according to the present invention.

FIG. 7 b is a flow chart showing a procedure for passing a strip.

FIG. 7 c is a flow diagram showing a block.

FIG. 7 d is a block diagram showing codebook processing.

FIG. 7 e is a diagram showing block delineation within a picture.

FIG. 7 f is a diagram showing codebook selection by content region.

FIG. 7 g is a diagram showing strip delineation according to region.

FIG. 7 h is a video sequence comprising a group of dependently codedpictures.

FIG. 7 i is a representation of the data composition of a scene.

FIG. 7 j shows syntax and semantic definitions of data elementsaccording to the present invention.

FIG. 8 a shows a conventional packetized elementary stream.

FIG. 8 b shows a private stream type within a multiplex.

FIG. 8 c shows conventional scenes and groups of pictures.

FIG. 8 d shows a conventional relationship coded frame and display frametimes.

FIG. 8 e shows codebook application periods.

DETAILED DESCRIPTION Overview of Tools

The PHD decoding process depicted in FIG. 4 b has two fundamental stagesof modem video coding. A first prediction phase 4130,1130 forms a firstbest estimate 4132,1135 of the target picture 4152,1175, using only theoutput state 4115,1115 of a base layer decoder 4110,1110 (and someminimal directives 4122), followed by a prediction error phasecomprising classification 4140,1120, enhancement decode 4120,1150 andapplication 4150 of correction 1165 terms that improve the estimate.

The overall PHD enhancement scheme fits within the template of theclassic spatial scalable enhancement architecture (FIG. 4 a). Therespective base layer decoders 4020,4110 are principally the same. Bothfundamental enhancement phases may operate concurrently in the receiver,and their respective output 4126,4032 added together at a later, thirdphase 4150, where the combined signal 4152 is sent to display, andoptionally stored 4160 for future reference 4172 in a frame buffer 4172.In a simplified embodiment the enhanced reconstruction 4152 may be sentdirectly to display 4162 to minimize memory storage and latency.

As part of the estimation phase 4130, the decoded base layer picture4115 is first interpolated according to parameters 4122 to match theresolution of the reconstructed HD image 4152. The interpolated image isa good first estimate of the target frame 4152. Traditionalinterpolation filters are applied in the preferred embodiment during theinterpolation process.

A first stage of the prediction error is to extract 4×4 blocks 1115 fromthe decoded base layer picture (4115) for classification analysis 4140.In order to keep computational complexity to a minimum, the preferredembodiment does not classify the interpolated base layer picture 4132,since the interpolated image nominally has four times the number pixelsas the base layer image 4115. The interpolated image 4132 is simply anenlarged version of the base layer image 4115, and inherently containsno additional information over the non-interpolated base layer image4115.

The preferred embodiment employs vector quantization to generatecorrection terms, in the form of 8×8 blocks 4126. Each block, orcodevector, within the codebook represents a small difference areabetween the interpolated predicted base image 4132 and the desired image4152. The codebook comprising VQ difference blocks are stored in a lookup table (LUT) 1160. The difference blocks are ideally re-used manytimes during the lifetime of the codebook.

Encoder

FIG. 5 c denotes the time order of the multi-pass base 5220 andenhancement layer (5230, 5240) video encoding processes. Nominally, thebase layer signal 5022 is first generated for at least the, period thatcorresponds to the enhancement signal period coded in 5230. Alternativeembodiments may jointly encode the base and enhancement layers, thusdifferent orders, including concurrent order, between 5210 and 5230 arepossible. The overall enhancement process has two stages: look-ahead5230 (FIG. 5 d) and real-time processes 5240 (exploded in FIG. 5 a). Theenhancement look-ahead period is nominally one scene, or access unitinterval for which the codebook is generated and aligned. The iterationperiod may be one scene, GOP, access unit, approximate time intervalsuch as five minutes, or entire program such as the full length of amovie. Only during the final iteration are the video bitstreams (5022,5252) actually generated, multiplexed into the program stream 5262, andrecorded onto DVD medium 5790. For similar optimization reasons, thefinal enhancement signal 5252 may also undergo several iterations. Themulti-pass base layer encoding iterations offer an opportunity in whichthe PHD look-ahead process can operate without adding further delays orencoding passes over the existing passes of prior art DVD authoring.

FIG. 5 b lists the databases maintained by the encoder 5110 look-aheadstages of FIG. 5 c. The enhancement codebook 5342 (database 5140) isconstructed by 5340 (described later) from training on blocks extractedfrom difference signal 5037 (database 5130). The codebook is lateremitted 5232, packed 5250 with other enhancement sub-streams (5234,5252) and data elements and finally multiplexed 5260 into the programstream 5262. In the preferred embodiment, the difference signal 5037 isgenerated just-in-time, on a block basis, from delayed pre-processedsignal 5010 stored in buffer 5013 (database 5160). Likewise, the baselayer signal 5032 (database 5120) is generated just in time from decodedSD frames (database 5150). Alternative embodiments may generate anycombination of the signals that contribute to the enhancement streamencoding process, either in advance (delayed until needed by buffers),or just-in-time.

The first two pre-classification stages 5310, 5320, described later inthis document, produce two side information arrays (or enhancementstreams) 5325 and 5315 (database 5180) that arc later multiplexed, alongwith the codebook, into the packed enhancement stream 5252. The resultsof the third pre-classification stage 5332 of FIG. 5 d may betemporarily maintained in encoder system memory, but are used only forcodebook training.

Although original HD frames (signal 5007) are in the preferredembodiment are passed only to the pre-processor 5010, furtherembodiments may keep the frames (database 5170) for multi-pass analysisin the classification or codebook training phases.

Run-time operations 5240, whose stages are detailed in FIG. 5 a, can begenerally categorized as those enhancement stages that produce packedbitstream elements for each coded enhancement picture. The enhancementdata may be buffered 5820 or generated as the final DVD program streamis written to storage medium 5790 master file. Buffering 5820 allows theenhancement stream to have variable delays to prevent overflow in thesystem stream multiplexer 5260. Enhancement may be generated in stepwith the base layer 5020 encoder at granularities of a blocks,macroblocks, macroblock rows and slices, pictures, group of pictures,sequences, scenes or access units. An alternate embodiment (FIG. 5 j) isto generate the enhancement stream 5252 after the base layer signal 5022has been created for the entire program, as would be the case if theenhancement is added to a pre-existing DVD title.

A second alternate embodiment is to generate the base and enhancementlayers jointly. A multi-pass DVD authoring strategy would entail severaliterations of each enhancement look-ahead process, while the joint baseand enhancement rate controllers attempt to optimize base andenhancement layer quality.

For best coding efficiency, the applied codebook and enhancement streamare generated after the scene, GOP (Group of Pictures), or otherinterval of access unit has been encoded for the base layer. The delaybetween base layer and enhancement layer steps is realized by buffers5013 and 5023.

The pre-processor 5010 first filters the original high-definition signal5007 to eliminate information which exceeds the desired rendering limitof the PHD enhancement process, or patterns which are difficult torepresent with PHD. The outcome 5012 of the pre-processor represents thedesirable quality target of the end PHD process. Film grain and otherartifacts of the HD source signal 5007 are removed at this stage.

The SD source signal 5017 is derived from the pre-processed HD signal5012 by a format conversion stage 5015 comprising low-pass filters anddecimators. The SD signal 5017 serves as source input for MPEG-2encoding 5020.

MPEG-2 encoder 5020 produces bitstream 5022, that after delay 5023, ismultiplexed as a separate elementary stream 5024 in the program streammultiplexer 5280.

The SD signal 5027 reconstructed by MPEG-2 decoder 5025 from delayedencoded SD bitstream 5024 is interpolated 5030 to serve as theprediction for the target HD signal 5014.

The prediction engine 5030 may also employ previously enhanced frames5072 to form a better estimate 5032, but nominally scales each picturefrom SD to HD dimensions.

The difference signal 5037 derived from the subtraction 5035 of thepredicted signal 5032 from the HD target signal 5014 serves as both atraining signal and enhancement source signal for the PHD encodingprocess 5050. Both source signals require the corresponding signalcomponents generation within the PHD encode process 5050 and enhancementcoding

The classifier 5040 analyzes the decoded SD signal 5027 to select aclass 5047 for each signal portion, or block, to be enhanced by the PHDencoding process 5050. The encoded enhancement signal 5052 is decoded bythe PHD decoder 5060, which in the encoder system can be realized as alook up table alone (5061) since the indices exist in pre-VLC (VariableLength Coding) encoded form within the encoder. The decoded enhancementsignal 5062 is added by 5065 to the predicted HD signal 5032 to producethe reconstructed HD signal 5067. The goal of the PHD encoder is toachieve a reconstruction 5067 that is close to the quality of the targetHD signal 5014.

The reconstructed HD signal 5067 may be stored and delayed in a framebuffer 5070 to assist the interpolation stage 5030.

The encoded PHD enhancement signal 5052 is multiplexed 5260 within theDVD program stream as an elementary stream with the base layer videoelementary stream 5024.

Some stages of the run-time operations are common to both the encoderand decoder. The encoder explicitly models decoder behavior when adecoded signal is recycled to serve as a basis for prediction 5072 infuture signals, or when the decoder performs some estimation work 5040of its own. For similar reasons, the MPEG-2 encoder 5020 models thebehavior of the MPEG-2 decoder 5025.

Pre-Processor (5010)

The primary responsibility of the pre-processor 5010 is to performformat conversion that maps the master source signal 5007 to the samplelattice of the HD target signal 5014.

The most common source format for HD authoring is SMPTE 274M, with 1920luminance samples per line, and 1080 active lines per frame. In order tomaintain a simple 2:1 relationship between the base and enhancementlayers, and to set a realistic enhancement target, the preferredenhancement HD coding lattice is twice the horizontal and verticaldimensions of the coded base layer lattice. For “NTSC” DVD's, this is1440×960 and 1408×960 for respective 720×480 and 704×480 base layerdimensions. For “PAL” DVD's with 576 active vertical lines, theenhancement dimensions are 1440×1152 and 1408×1152 respectively. Thebase layer will assumed to be 720×480 for purposes of this description,although the enhancement process is applicable to any base andenhancement dimension, and ratio.

A skilled engineer can choose from many image scaling designs, includingwell known poly-phase FIR filters, to convert the first 1920×1080 framelattice of 5012 to the second 1440×960 lattice of 5017. Another possibleformat for either or both of the input 5012 and output 5017 sides isSMPTE 296M, with 1280×960 image dimensions. A corresponding formatconversion stage 1482 in the decoder maps the PHD coded dimensions tothe separate requirements of the display device connected to displaysignal 1482. Common display formats include SMPTE 274M (1920×1080×30i)and SMPTE 296M (1280×720×60p).

General format conversion pre-processing essentially places the targetsignal in the proper framework for enhancement coding. The goal ofpre-processing is to produce a signal that can be efficientlyrepresented by the enhancement coding process, and assists theenhancement coder to distribute bits on more visibly important areas ofthe picture. Several filters are employed for the multiple goals ofpre-processing.

A band-pass filter eliminates spatial frequencies exceeding a user orautomatically derived target content detail level. The band-pass filtercan be integrated with the format conversion scaling filters. The formatscaling algorithm reduces the 1920×1080 HD master format to the 1440×960coding format, but additional band-pass filtering smooths the contentdetail to effectively lower resolutions, for example, 1000×700.

Adaptive filtering eliminates patterns that are visually insignificant,yet would incur a bit cost in latter encoding stages if left unmodifiedby the pre-processor. Patterns include film grain; film specs such asdirt, hair, lint, dust;

A classic pattern and most common impediment to efficient coding issignal noise. Removal of noise will generally produce a cleaner picture,with a lower coded bit rate. For the PHD enhancement process, noiseremoval will reduce instances of codebook vectors that would otherwisebe wasted on signal components chiefly differentiated by noise. Typicalnoise filters include 2D median, and temporal motion compensated IIR andFIR filters.

Downsample (5015)

The base layer bitstream complies with MPEG-2 Main Profile@Main Levelvideo sequence size parameters fixed by the DVD specification. AlthoughMPEG-2 Main Profile@Main Level can prescribe an unlimited number ofimage size combinations, the DVD specification limits the MPEG-2 codingparameters to four sizes (720×480, 704×480, 720×576, and 704×576), amongwhich the DVD author can select. The DVD MPEG-1 formats (352×240 and352×288) are not described here, but are applicable to the invention.The HD target sample lattice 5012 is decimated 5015 to the operationallattice 5017 of the MPEG-2 5020. Downsampling 5015 may be bypassed ifthe encoder 5020 is able to operate directly upon HD formats, forexample, and is able to perform any necessary conversion to the DVD baselayer video format. In prior art, downsampling 5015 will execute masterformat conversion, such 24p HD (SMPTE RP 211-2000) to the SD formatencoded by 5020.

Downsampling may be performed with a number of decimation algorithms. Amulti-tap polyphase FIR filter is a choice.

MPEG-2 Encoder (5020)

The MPEG-2 encoder 5020 nominally performs as prior art encoders for DVDauthoring. Although the invention can work with no changes to the baselayer encoder 5020, improvements to the overall reconstructedenhancement layer video can be realized through some modification of thebase layer encoding process. In general, any operation in the base layerthat can be manipulated to improve quality or efficiency in theenhancement layer is susceptible to coordination with the enhancementprocess. In particular, operation of the DCT coefficient quantizermechanisms quant_code and quantization_weighting_matrix .can becontrolled to maintain consistent enhanced picture quality. In somecombinations of base and enhancement data, this would be more efficientthan applying additional bits to the corresponding area in theenhancement layer. In an advanced design, the rate control stage of theencoder 5020 could have dual base and enhancement layer rate-distortionoptimization.

Improved motion vectors coding in the base layer may benefit modes ofthe enhanced prediction stage 5030 that employ motion vectors extractedfrom the base layer signal 5022 to produce interpolated predicted frames(a feature of an alternate embodiment described later in thisspecification). Motion vector construction is directly operated byrate-distortion optimization with feedback from both the base andenhancement reconstruction.

The encoder may also need to throttle back the bitrate to ensure thecombination of enhance and base bitstreams do not exceed DVD buffercapacity.

Prediction (5030)

The prediction scheme forms a best estimate of the target signal bymaximizing use of previously decoded data, and thereby minimizing theamount of information needed for signaling prediction error. For theapplication of picture resolution and detail enhancement, a goodpredictor is the set of image interpolation algorithms used in scalingpictures from one resolution, such as an intermediate or coded format,to a higher resolution display format. These scaling algorithms aredesigned to provide a plausible approximation of signal content sampledat higher resolution given the limited information available in thesource lower resolution picture.

Overall, the base layer decoded image 6110 extracted from signal 5027 isscaled by a ratio of 2:1 from input dimensions 720×480 to an outputdimension of 1440×960 of the signal 5032 to match the lattice of thetarget 5014 and enhanced images 5067 so that the predicted signal 5032image 6120 may be directly subtracted 5035 from the target signal 5014,and directly added 5065, 6130 to the enhancement difference signal 5062image 6140 to produce the enhanced picture 6150. Other ratios and imagesizes are applicable. In some picture areas or blocks, the predictedsignal 5032 is sufficient in quality to the target signal 5014 that noadditional information 5052 need be coded.

The order of the stages within the prediction 5030 function of thepreferred embodiment is depicted in FIG. 6 a. Other orders are possible,but the preferred order is chosen as a balance between implementationcomplexity and performance, and for dependencies with the base layerbitstream such as the de-blocking stage's use of quantizer step sizes.Starting with the base frame 6010, 6110 extracted from signal 5027, ade-blocking filter 6020 is applied to reduce coding artifacts present inthe base layer. Although good coding generally yields few artifacts,they may become more visible or amplified as a result of the scalingprocess 6030, or plainly more visible on a higher definition screen.De-blocking reduces unwanted patterns sometimes unavoidably introducedby the MPEG-2 base layer encoding process 5020.

The de-blocking filter of ITU-T H.263 Annex J is adapted to 6020. Somestages of the Annex J filter require modifications in order to fit theinvention. For example, the de-blocking filter is performed as apost-processing stage after the image has been decoded, not as part ofthe motion compensated reconstruction loop of the base layer decoder.The quantization step function is remapped from the H.263 to the stepsof the MPEG-2 quantizer. The strength of the de-blocking filter isfurther regulated by a global control parameter transmitted with eachenhanced PHD picture. The PHD encoder sets the global parameter toweight the Annex J STRENGTH constant according to analysis of thedecoded picture quality. Since the quantizer scale factor is not alwaysan indication of picture quality or coding artifacts, the PHD encoderaims to use the global parameter to set the STRENGTH value to minimalfor pictures with excellent quality, thus de-blocking is effectivelyturned off when it is not needed or would do unnecessary alterations tothe picture.

A poly-phase cubic interpolation filter 6030 derives a 1440×960 image6035 from the de-blocked standard definition 720×480 image 6025.

In an alternative embodiment (FIG. 6 c functional blocks and FIG. 6 ddata blocks), a subset of pictures within a sequence or GOP arealternatively predicted from a combination of previously decoded baselayer and enhanced pictures 6320,6322 stored in frame buffer 6225—asubset of frame buffer 5070. This variation of a predicted enhancementpicture is henceforth referred to as a temporally predicted enhancementpicture (TPEP) 6345. TPEP resembles the B-frame or “bi-directionally”predicted frames since they borrow information from previously decodedframes that in display order are both future and past. The differenceenhancement 6320, 6322 from previously decoded pictures is reapplied tothe current picture 6315 as a good estimate of the enhancementdifference 6140 that would be otherwise transmitted as enhancement datain non-TPEP pictures. TPEP is a tool for reducing the overall or averagebitrate of the enhancement layer since data is not often coded for TPEPblocks. If difference mode is enabled in the header of TPEP pictures, a1-bit flag prefixes each TPEP block indicating whether differenceinformation will be transmitted for the block. TPEP pictures are enabledwhen the corresponding base layer picture is a B picture; the scaledmotion information 6235 from the base layer picture instincts the MCP6235 to create the prediction surface 6325 that is combined 6340 withthe interpolated base frame 6315.

Classification

While Standard Definition (SD) and High Definition (HD) images capturedof the same scene differ superficially by the density and size of theirrespective sample lattices (1440×960 vs. 720×480), they maysubstantively differ in content, in particular when analyzed in thefrequency domain. Generally, a hierarchical relationship should exist inthat the information in the SD image is a subset of the HD image, suchthat the SD image may be derived from the HD image through operationssuch as filtering and sub-sampling. (Eq.1)

SD=sub-sample(HD)  (Eq. 1)

In the spatial domain, an HD image can be represented as the sum of afirst base image (B) and a second difference (D) image:

B=sub-sample(HD)

D=HD−B  (Eq. 2)

HD′=B′+D  (Eq. 3)

In this example, the difference image (D) contains the high frequencycomponents that distinguish the HD image from the SD image, while thebase image (B) contains the remaining low frequency information. Whenthe base image (B) by itself can serve as the SD image, the differenceimage (D) could then be formulated to contain the set of informationthat is present only in the HD image, not the SD image.

Further, the SD image can be sampled at a reduced resolution, with asmaller lattice (such as 720×480), sufficient to contain the lower halfof the frequency spectrum, and later scaled (SD′) to match the samplelattice (e.g. 1440×960) of the HDTV image where it may be easilyrecombined in the spatial domain with the difference image (D) toproduce the reconstructed HD image (HD′).

While the lower frequencies are significantly more important than highfrequencies in terms of perceptible contribution to the overall image(HD′), the high frequency information is still needed to establish the“look and feel” of an HD image.

Although the difference image may be expected to contain up to threetimes more information than the base image, not all portions of thedifference image contribute equally to the overall perceptible qualityof the final reconstructed HD image. The essential information in (D)needed to emulate the look and feel of the HD image may in fact be asmall subset of D, in particular concentrated along edges and areas oftexture, and may be further approximated very coarsely. This concept isessentially supported by the practice in the block coding methods ofJPEG and MPEG where high frequency DCT coefficients are more coarselyquantized than low frequency DCT coefficients.

The MPEG coding tools are not optimized for coding these essentialdifference areas efficiently at extremely low bit-rates (or in otherwords, high compression factors). MPEG is tuned towards visualapproximation of an image with a balance of detail and generic contentat appropriately matched resolutions. For example, the luminance samplesof a typical still frame will be represented as an MPEG intra-frame (I)in approximately one fourth the rate of the “non-coded” PCM frame, andthe average predicted frame (P,B) only one fifteenth the size of the PCMframe.

The classifier stage of the invention serves as a key tool foridentifying those areas of the picture of greater subjective importance,so that enhancement coding may be emphasized there. At the same time,the process also objectively places emphasis on those areas where thedifference energy is greater, such as edges.

Strong horizontal, vertical, and diagonal edges, for example, can beidentified at lower resolutions, such as the SD base layer within the SDimage. It is possible to identify areas that should result in acombination of high frequency and high perceptible patterns in the HDimage. Unfortunately, sufficient clues in the base image are notaccessible to accurately estimate the actual difference information forthose areas, although reasonable guesses bounded by constraintsimprinted in the base layer are possible, and have been developed byvarious prior “sub-pixel” developments. To meet real-time implementationconstraints, prior art interpolation schemes would generate “synthetichighs” through contrast enhancement or sharpening filters. The mostcommon algorithm for interpolating image is a filter that convolves thelower resolution samples with a curve that models the distribution ofenergy in the higher resolution sample lattice, such as the sinc( )function.

Superficially sharp, high resolution images restored by syntheticallymeans from low resolution images often looks contrived or as anartificial byproduct, and quality gains may be inconsistent.

Accurate identification of picture areas is possible with knowledge ofthe original HD image, but such an image is available only to theencoder residing at the transmitter side. Enhancement information can beexplicitly transmitted with this knowledge to guide the HDreconstruction process, and thus produce more natural looking “highs”.However enhancement data can easily lead to a significant bit rateincrease over the base layer data.

The more accurate the highs can be estimated by the receiver, the lessenhancement information is needed to improve the reconstructed HD signalto a given quality level. A particular tool useful for minimizing thevolume of enhancement information is classification.

Classification can be used to partially predict the enhancement layerand/or prioritize those areas that need to be enhanced. Classificationalso permits different coding tools to be used on different classes ofpicture data. For example, in flat areas the SD to HD interpolationalgorithm may dither, while pixels determined to belong to an edge classmay benefit from directional filtering and enhancement data.

As appropriate for the overall enhancement technique, classification canbe accomplished in the frequency or spatial domains. A classifier isalso characterized by the granularity of the classified result (such ason a per pixel or block basis), and by the window of support for eachgranule.

The window of the classifier is the size of the support area used in theclassification analysis. For example, to determine the class of a singletarget pixel, the surrounding 5×5 area may be measured along with thetarget pixel in order to accurately measure its gradient.

Familiar to video compression, a good balance between implementationcomplexity, bitrate, and quality can be achieved with block-basedcoding. The negative tradeoff is manifested by inaccuracies that resultat block edges and the other blocking artifacts.

The preferred PHD classification scheme employs block-based frequencyand spatial domain operators at a granularity of 4×4 pixels with respectto the base layer, and 8×8 pixels with respect to the HD image. Localimage geometry (flat, edge, etc.) is first determined through a seriesof comparisons of measurements derived from frequency coefficients of a4×4 DCT taken on a non-overlapping block within in the base image.Overlapping is also possible, but not implemented in the preferredembodiment. The small 4×4 block size has many of the desired propertiesof a local spatial domain operation, but with greater regularity andreduced complexity compared to both per-pixel granular operations, andgenerally most known effective all-spatial domain operations.

Calculating Classification Components

FIGS. 3 b and 3 d provide the fundamental stages of the preferredclassifier embodiment that are common to both the encoder and decoder.FIG. 3 d discloses the classifier' component calculations 3130 of FIG. 3b.

Blocking

Blocks of data are extracted from the input frame 3100 in the processingorder of the enhancement decoder. The preferred processor order israster, from left to right and top to bottom of the picture, withnon-overlapping blocks. Alternate embodiments may overlap blocks inorder to improve classification accuracy. For example, a 3×3 targetblock may be processed from a 4×4 input block. In the 3×3 within 4×4block example, the overlap areas would comprise a single row and columnof extra pixels. Each successive 3×3 picture area would then beprocessed from a 4×4 block with a unique combination of samples formedfrom the base picture. The 4×4 input block would step three pixels foreach advance in either or both the x and y directions. A new set ofclassification parameters would be derived for each 3×3 picture area.Other overlaps are possible, but in general, the overlap and targetblocks may be arbitrarily shaped as long as the base and enhancementlayers are aligned.

DCT

In the preferred embodiment, the DCT-II algorithm is applied in the 4×4DCT 3312 to produce the coefficients 3314 whose combinations are used asfeature component measurements 3332 for the decision stage 3140.Variations include the DCT-I and DCTIII, non-DCT algorithms, andpseudo-DCT algorithms such as those experimented with by the ITU-T H.264study group. Generally, any transform which produces coefficients usefulin the classification of a picture area can substitute for the preferredblock DCT, however adjustments to the ratio calculations in 3130 anddecision tree 3140 may be necessary to account for the differentcharacteristics of each transform's unique coefficient sets.

The 8-bit precision of the transform coefficients and 16-bitintermediate pipeline stages are sufficient to support the expansion ofdata in the transform size and the accuracy needed to discriminate oneclass from another. The preferred transform is designed to operatewithin the 16-bit SIMD arithmetic limitations of the Intel MMXarchitecture which serves as an exemplary platform for PHD DVDauthoring.

Spatial Analysis

The Weber function provides a more accurate measurement of picture areaflatness than a single combination of DCT coefficients.

The Weber component 3322 calculated in 3320 follows the formulasummarized as:

-   -   compute difference between max value of block and average block        value if the difference/average<=0.03, then it is flat        (isFlag=1), else isFlag=0.

Frequency Analysis

Component generator 3330 takes measurements 3132 conducted on the 4×4blocks and produces decision variables 3332, 3132 used in the decisionprocess 3140 to create classification terms 3142. The block measurements3132 comprise both frequency measurements 3314 (in the preferredembodiment realized by the 4×4 DCT transform 3312) and spatial domainmeasurements 3322 (in the preferred embodiment realized by a flatnessoperator 3320).

Input blocks 3310, 3122 formatted from the base layer reconstructedimage 3100 are transformed via the 4×4 DCT 3312, producing coefficients3314. The component generator stage 3332 takes sets of coefficients 3314shown in FIG. 3 e, and squares and sums coefficients within each set toproduce class components 3332, P1 through P7. Each set of DCTcoefficients, and its resulting measurement term (P1 . . . P7),represents the identifying characteristic of a geometric shape such asan edge, texture, flat area.

The seven 4×4 DCT coefficient templates in FIG. 3 e shows increasinghorizontal frequency is along the U-axis with set of indices {0, 1, 2,3}, and increasing vertical frequency along the V-axis with indices{A,B,C,D}.

Each of the components P1 . . . P7 represent the following geometryfeatures: P1—horizontal edges, P2—horizontal texture, P3—vertical edges,P4—vertical texture, P5—diagonal edges, P6—texture, andP7—energy/variance of the block.

(P1) diag=Bl*B1+C2*C2+D3*D3

(P2) inf0=BO*B0+CO*C0+DO*DO+C1*C1+D1*D1+D2*D2

(P3) infl BO*B0+CO*C0−+DO*DO

(P4) sup( )Al*A1+A2*A2+A3*A3+B2*B2+B3*B3+C3*C3

(P5) supi−Al*A1+A2*A2+A3*A3

(P6) text−C2*C2+C3*C3+D2*D2+D3*D3

(P7) tot=diag+sup( )+into

Ratios

From the seven component measures (P1 . . . P7), eight ratios (R0 . . .R7) arc derived that are used in the decision process 3140 to select theclass for each block.

RO=diag/tot

R1=sup( )/(sup°+infO)

R2=supl/sup°

R3=info/(sup( )+info)

R4=infl/inf 0

RS=text/(sup°+infO)

R6=supl/(sup°+infO)

R7=infl/(sup°+infO)

Pre-Calculated Ranges

In order to improve accuracy of the codebook and run-time classificationpasses, two pre-classification passes 5310, 5320, 5330 are made throughthe decoded base layer signal 5027, 5305, to measure the statistics ofclassification components. Specifically, thresholds 5317 and energyranges 5327 are produced in the first and second passes respectively.The third classification pass 5330 selects the class for each trainingblock 5332 used in codebook generation stage 5340. The codebook istrained on the decoded base layer signal; the results of the thirdpre-classification stage therefore 5332 model (sans IDCT drift error)the run-time classifier 5040 results of downstream decoder classifier.

Ratios RO . . . R7 are calculated in the classification stage as above,and then compared to pre-determined thresholds to establish 17 energyranges 5327.

Ranges and thresholds (shown collectively as side information 5234) aremaintained in memory 5180 for later application in the class decisionstage 3140. To save computation time, and spare the decoder from havingto add significant latency, the encoder packs the ranges and thresholdsinto the PHD stream 5252, where on the receiver side, they are laterparsed and integrated into the state machine 3620 by the PHD decoderduring each codebook update.

To improve accuracy of classification, the components used in theclassification decision process are adaptively quantized accordingtraining block statistics. The quantized levels are indicated bythresholds 5315 which are calculated from an equi-probable partitioningof histograms measured during the first pre-classification training pass5310.

Pass 1, generate adaptive quantization thresholds: For each trainingblock.

if ((R1 > 0.60) && (R2 <= 0.90))   hist_add( hist1, R1 ); else if ((R1 >0.60) && (R2 > 0.90))   hist_add ( hist2, R1 ); else if ((R3 >0.60) && (R4 <= 0.90))   hist_add( hist3, R3 ); else if ((R3 >0.60) && (R4 > 0.90))   hist_add( hist4, R3 );

Hist_add(argl, arg2) updates respective histogram (indicated by argl)with the data point arg2. Each histogram is allocated to track a rangeof values divided into a specified number of partitions. Each update ofarg2 will increment the corresponding partition identified by arg2 byone count.

At the end of the training sequence, hist_convg(argl, arg2, arg3, arg4)partitions thresholds 5315 (arg3) into arg4 number of equi-probablepartitions according to the statistics stored in the respectivehistogram argl

At the end of the training session.

-   -   hist_convg(histl, hcenters, threshl, 2);    -   hist_convg(hist2, hcenters, thresh2, 5);    -   hist_convg(hist3, hcenters, thresh3, 2);    -   histconvg(hist4, hcenters, thresh4, 5);        The second parameter, arg2, of Hist_conv( ) provides additional        statistics including the average and standard deviation squared        of each partition.

Pass 2, Measure Energy:

Note: is Flat is the result of the Weber calculation 3320.

if (isFlat)  idx = 0; else {  if (RO >= 0.55)   idx = 1;  else  {   if((R1 > 0.60) && (R2 <= 0.90))   {    if (R1 < threshl [ 0] )     idx =2;   else     idx = 3;   }   else if ((R1 > 0.60) && (R2 > 0.90))   {   if (R1 < thresh2 [0] )     idx = 4;    else if (R1 < thresh2 [1] )    idx = 5;    else if (R1 < thresh2 [2] )     idx = 6;    else if (R1< thresh2 [3] )     idx = 7;    else     idx = 8;    }   else if ((R3 >0.60) && (R4 <= 0.90))   {    if (R3 < thresh3[0])     idx = 9;    else    idx = 10;   }   else if ((R3 > 0.60) && (R4 > 0.90))   {     if(R3 <thresh4[0])      idx = 11;     else if (R3 < thresh4[1])      idx = 12;    else if (R3 < thresh4[2])      idx= 13;     else if (R3 <thresh4[3])      idx = 14;     else      idx = 15;   }   else     idx =16; t [idx] [count [idx] ] = Etot ; count [idx] = count [idx] + 1;min_energy_class [idx] =  MYMIN ( min_energy_class [idx] , Etot ) ;max_energy_class [idx] =  MYMAX ( max_energy_class [idx] Etot );

At the end of the second pre-classification pass 5320 of the trainingsequence, the statistics in temporary variable arrays t[ ] and count[ ]are used to calculate 17 energy_range[ ] 5325 constants used in theclassification stage.

for (i = 0; i < 17; i++) {   median (count [i] , &t [i] [0] ,&median_val) ;   energy_(—) range [i] = median_val ; }

Determining Class by Decision Tree

To arrive at a specific class, the classifier uses the componentmeasurements produced in 3510, 3330, to descend a decision tree,comparing class components 3332 and pre-calculated ranges (3102, 5180,5240, 5234, 5315). The generic cyclical flow of the classificationprocess is given in FIG. 3 f. Comparisons are made 3520 until a stateprocess indicates that a class has been arrived at 3530. With the binarydecision branch process depicted, the number of iterations should beapproximately the logarithm of the number of available classes. Means ofimplementing the decision tree include procedural code (nested ifstatements) given below, and parallel flow-graph testing (not shown).

A state machine realization of the decision tree is given in flowchartFIG. 3 g. The state machine is expected to be the easiest Stateparameters table 3620 is indexed by variable L, initialized to zero3610. The resulting state parameters 3621 include branch positiveaddress L1, branch negative address L2, classification componentidentifiers pl and p2, multiplier constant k, offset T, and terminationbits el and e2.

Component identifiers pl and p2 select which classification ratios inthe set Pl . . . P7 are to be compared in 3640. The values for p1 and p2are selected 3630 from the class component register array cc andcompared as a and b in formula 3640. The branch addresses L1 are thenext location in the state code 3620 that the state program reaches ifthe comparison in 3640 is positive, and L2 is the location if thecomparison is negative. If either or both of the comparison resultsindicate a terminal condition, that is a terminal node with a specificclass is finally reached, then either or both terminal state bits el, e2will be set to ‘1’ potentially causing the loop to exit Y at 3650. Inaterminal cases (where E==1), state variables L1 and L2 encode the classindex 3632 which forms part of the state 3142 in FIG. 3 b needed toperform, at 1 east, the LUT 3150.

A procedural example of the decision tree is below. Energy_class:

if (isFlat)  energy_class[i] = 0; else {  if (RO >= 0.55) // diagonal  {  if (Etot < energy_range[1])     {       energy_class[i] 1;     }  else     {       energyclass[i] = 2;     }  }  else   {     if ((R1 >0.60) && (R2 <= 0.90))       {         if (R1 < thresh1[0]) //vert_text_0         {           if (Etot < energy_range[2])            energy_class[i] 3;           else            energy_class[i] = 4;         }         else // vert_text_1        {         if (Etot < energy_range[3]) // vert_text          energy_class[i] = 5;         else           energy_class[i] =6;       }   }   else if ((R1 > 0.60) && (R2 > 0.90))   {   if (R1 <thresh2[0]) // count_vert_0     {       if (Etot < energy_range[4])        energy_class[i] = 7;       else         energy_class[i] = 8;    }   }   else if (R1 < thresh2[1]) // vert_1   {     if (Etot <energy_range[5])         energy_class[i] = 9;       else        energy_class[i] = 10;   }   else if (R1 < thresh2[2]) // vert2  {     if (Etot < energy_range[6])         energy_class[i] = 11;      else         energy_class[i] = 12;   }   else if (R1 < thresh2[3])// vert_3   {     if (Etot < energy range[7])         energy_class[i] =13;       else         energy_class[i] = 14;   }   else // vert_4   {    if (Etot < energy_range[8])         energy_class[i] = 15;       else        energy class[i] = 16;   }   else if ((R3 > 0.60) && (R4 <=0.90))   {     if (R3 < thresh3[0]) // text_0     {       if (Etot <energy_range[9])           energy_class[i] = 17;         else          energy_class[i] = 18;     }     else // horz_text 1     {    if (Etot < energy_range[10])         energy_class[i] = 19;      else       energy_class[i] = 20;     } } else if ((R3 0.60) &&(R4 > 0.90)) {   if (R3 < thresh4[0]) // horz_0   {     if (Etot <energy_range[11])       energy class[i] = 21;     else      energy_class[i] = 22;   }   else if (R3 < thresh4[1]) // horz_1  {   if (Etot < energy_range[12])       energy class[i] = 23;     else      energy_class[i] = 24;   }   else if (R3 < thresh4[2]) // horz_2 {if (Etot < energy_range[13])     energy_class[i] 25;   else    energy_class[i] = 26; } else if (R3 < thresh4[3]) // horz_3 { if(Etot < energy_range[14])   energy_class[i] = 27; else   energy_class[i]28; } Else {   if (Etot < energy_range[15]) // horz_4      energy_class[i] = 29;     else       energy_class[i] = 30;        count ++;     }   }   else // ((R5 < 0.35) && (R6 < 0.65) && (R7< 0.65))   { // text_0     if (Etot < energy_range[16])        energy_class[i] = 31;       else         energy_class[i] = 32; }

Entire scenes, or individual pictures often do not contain significantdetail in the original high-definition format signal beyond the detailthat would be prescribed in any standard definition derivative of thehigh-definition signal. In such cases when there is insufficientdifference between the high definition original signal 5012 andpredictive signal 5032, it is more efficient to turn off enhancementblock coding, while predictive interpolation continues to operate underboth conditions in one mode or another.

To determine whether enhancement blocks should be sent for an area(encapsulated as a stripe), picture, or scene, the selective enhancementanalyzer 5420 estimates the perceptivity of the difference signal 5037for each block prior to both the VQ codebook training and run-timecoding phases. Although many models exist for perceptivity, the simpleenergy formula calculated as the square of all N elements within theblock serves as a reasonable approximation. The preferred embodimentapplies the following formula:

$e = {\sum\limits_{i = 0}^{N - 1}\left( {{block}\lbrack i\rbrack} \right)^{2}}$

Three control parameters 5422 regulate the selection algorithm in 5420.The first user control parameter, energy_threshold, sets the level ofenergy for a block to meet in order to be selected for enhancement bythe encoder. Since the measurement is made on the difference signal5037, only the encoder can make such a judgment, although special casessuch as flat areas (described earlier) that do not have associatedindices are determined by the receiver through measurements on the baselayer signal.

User control parameter stripe_block_ratio_threshold sets the minimumratio of selected blocks within a stripe that must meet the perceptivitycriteria in order for the slice to be coded. User control parameter block_max sets the level in which, regardless of the ratio of selectedenhancement blocks, the stripe would be coded. This accounts forisolated but visually significant blocks.

Stripe headers include a 3-bit modulo index strip_count er so that thedecoder can distinguish between non-coded gaps in the enhancementpicture and stripes that have been lost to channel loss such as droppedor corrupted packets.

Blocks that do not meet the enhancement threshold are not applied duringthe VQ training process.

The is_picture_enhanced variable in the picture header signals whetherenhancement blocks are present for the current picture. For fingergranular control, the is_strip_enhanced flag in the strip header canturn enhancement blocks on or off for all blocks within a strip( ). Inmany cases, only a small subset of the picture has sufficient detail tomerit enhancement, usually those areas that the camera had in focus. Insuch cases, the encoder can adapt the strip( ) structure to encapsulateonly those detail areas, and leave the rest of the picture withoutstrip( ) coverage. The x-y position indicators within the strip( )header allow the strip( ) to be positioned anywhere within the picture.

PHD Run-Time Encoding (5050)

Enhancement data 5052 is generated for those blocks whose class hasassociated enhancement blocks 5062. Of the thirty three classes, class0, the category for flat areas, requires no transmission of indices. Thestatistical expectation is that at least one in three blocks will beclassified as flat, and for some scenes, flat blocks will constitute amajority of blocks. Thus the bitrate savings can be substantial by nottransmitting enhancement block indices for areas that do notsufficiently benefit from enhancement. Since the encoder and decoderhave an identical understanding of the enhancement syntax and semantics,the decoder parser does not expect indices for non-coded enhancementblocks.

For those classes with associated enhancement data, the VLC index ispacked within the enhancement bitstream 5262 along with otherenhancement elements. The combination of class and the VLC index arc allthat is needed to perform an enhancement pattern lookup 5060, where adifference block is generated 5062 and added 5065 to the correspondingpredicted-interpolated block 5032. The same lookup procedure isperformed in the receiver.

Small discrepancies in the reconstructed enhanced signal 5067 may existdue to difference among standard-compliant MPEG video reconstructions5024. No one model of the decoder 5025 applies universally. Drift freereconstruction is possible only if the IDCT in the encoder is matched tothe IDCT in the receiver. The difference signal, or drift, between themodel decoder 5025 and the actual downstream decoder originates due toround-off errors in the integer approximation of the standard-definedfloating point IDCT algorithm. The drift should be limited to anoccasional least significant bit difference, but in pathological casesdesigned to accumulate worst case patterns, drift has been known tobuild visible artifacts. Consequentially, drift can cause discrepanciesbetween the encoder model classifier result 5047 and classificationresult 4142 in the downstream decoder. With proper threshold design,these discrepancy cases are rare and detectable through. theclass_checksum mechanism in the header of each strip( ). Whenclass_checksum and the receiver calculated checksum differ, enhancementis not applied for those blocks for which the checksum applies. Thespecific class_checksum element applies to all blocks contained withinthe strip°.

The preferred embodiment applies the well known CRC-32 algorithm togenerate the bitstream checksum class_checksum and receiver checksum towhich it is compared. Other hash algorithms could be applied, but CRC-32circuitry is common in existing receivers with MPEG-2 video decoders.

Entropy Coding

The JPEG-2000 arithmetic coder is utilized by the invention for bothcodebook and enhancement block index transmission.

New codebooks are transmitted as raw samples. One codebook is sent foreach class that has specified transmitted indices. For classes that donot have codevectors, the size_of class variable (FIG. 7 j) is set tozero. The order of the codevectors within each codebook is at thediscretion of the encoder. The encoder should take care that the indicescorrespond to the correct codevector entry within the transmitted ordercodebook table.

-   -   Cbk[class_num][k=sample (8 bits);

Codebook updates arc sent as run-length encoded differences betweencorresponding blocks in the first codebook and the second codebook. Oneset of context models are created for each class. A first context modelmeasures run of zeros, while the second context addresses amplitude.

-   -   Diff_cbk [c] [v] [k]=new_cbk [c] [v] [k]−prev_cbk [c] [v] [k]

The difference codebook, cliff cbk, is calculated as the sample-wisedifference between the new codebook, new_vector, and the old codebook,prey cbk. Most cliff cbk samples will be zero, followed by smallamplitudes.

Specific arithmetic coding context models are created for each class ofthe enhancement block indices. The first context is the original indexalphabet to each class sub-table. A second context is the average of thepreviously transmitted above and left blocks.

The arithmetic coder is reset for each strip.

PHD Decoding

PHD decoding is a subset of the encoder operation, and is preciselymodeled by the encoder as illustrated in FIG. 5 a. Specifically, MPEG-2decode base layer 5025 is 4110, predictive interpolation 5030 is 4130,classifier 5040 is 4140, VQ decoder 5060 is 4107, adder 5065 is 4150,and frame buffer store 5070 is 4170.

Codebook Generation

Virtually any codebook design algorithm can be used to generate theenhancement codebook 5140. The codebook could also be selected from aset of universal codebooks rather than created from some trainingprocess on the video signal to be encoded. The preferred PHD vectorquantization codebook design algorithm is a hybrid of the GeneralizedLloyd Algorithm (GLA), Pair-wise Nearest Neighbor (PNN), and BFOSalgorithms described in [Garrido95]. The hybrid is continuously appliedto each video scene. Training sequences 5130 are derived from a set offiltered HD images 5160, 5012, rather than original HD images 5007,5170. Although it would be less expensive not to have the pre-processingstage 5010, the original HD source images are not used for comparisonsince it may contain data patterns that are either unnecessary for theapplication, or unrealistic to approximate with PHD coding. Thedifference signal 5332, 5037 generated as the difference between thecleaned signal 5014 stored in 5013, 5160 and the interpolative-predictedsignal 5032 is then fed to the codebook generator 5340.

A potential codebook 5140 is transmitted along with each scene, where itis then parsed by the PHD decoder at the receiver side, and stored inlong term memory 5160 for application throughout a scene or, in specialcases, applied repeatedly in future scenes. Syntax

The PHD syntax is structured to a hierarchy (FIG. 7 e) resemblingtraditional video coding layers known for efficient and robust parsing.A scene roughly corresponds to a typical video sequence (FIG. 7 h), andin addition to codebook updates, includes the energy thresholdparameters 5317, 5327 used in the classification stages. Picture headersenhancement_picture( ) delineate sets of indices corresponding to theenhancement blocks for a given picture. The picture header identifiesthe current enhancement picture number, picture_number, and the picturepayload comprises one or more strips that select which codebook codebooknumber is to be applied for those blocks contained within the strip.

Referencing Multiple Codebooks

Duration of Codebook:

A codebook is created for application upon a scene which typically lastsfrom half a second to several seconds, such as 8210, 8220, and 8230depicted in FIG. 8 c. In extreme cases, the lengths of scenes can rangefrom a few pictures to several minutes (thousands of pictures). Sinceevery scene has unique image statistics and characteristics, codebooksoptimized for each scene will produce better quality results for a givenindex rate. The overhead of sending codebooks also significantly impactsthe quality-rate tradeoff. Frequent transmission of codebooks willoffset the index quality gains, and potentially penalizing quality inthe base bitstream (if the base stream is jointly optimized), or leaveless room for future codebooks on the disc volume. Some scene changes,such as camera angle cuts with similar background (e.g. two characterstalking to each other) may precipitate codebooks that largely overlapwith previously sent codebooks. The differential and dynamic codebookupdate mechanisms disclosed herein address these cases. Pointers topreviously sent codebooks (FIG. 8 e) may also be more efficient forshort, repeating scenes.

The PHD advantage of exploiting long-term correlations is partlyillustrated in FIG. 8 c by the ability of a codebook (aligned to ascene) to span periods exceeding the nominal enforced “group ofpictures” (GOP) dependency periods, and thus saves bits compared to astrategy where a codebook is automatically sent for each GOP. Thus, forexample instead of transmitting a codebook every 0.5 seconds—the periodof the Intra-picture or GOP—the codebook need only be transmitted everyfew seconds. The random access period for the enhancement layer willthus consequently be greater than the base layer, but as long as a baselayer picture can be built with the normal short latency, a goodapproximation for the purposes of non-predetermined trick play can besatisfied. New codebooks are forced by the DVD authoring tools forpre-determined jumps within the DVD volume, such as angle or layerchanges. Thus playback along the pre-constructed linear stream timelinewill maintain constant enhanced picture quality.

In this invention, GOP is applied more widely to mean independentlydecodable collection of pictures, typically constructed in MPEG videostream to facilitate random access and minimize DCT drift error.“group_of_pictures_(—)0” has a narrower meaning in the MPEG videospecification than this description, but fits within the definitiongiven here. For this invention, GOP is a generic term, and superset ofthe formal MPEG definition, that delineates any collection ofdependently coded pictures. The duration of the GOP is typically 0.5seconds in DVD applications, but the exact boundary of a GOP may headjusted for scene changes or coding efficiency.

Random access to a codebook can be optimized for scene changes, bufferstatistics, chapter marks, and physical models such as location of thescene data within the disc.

Nominally, multiple bitstream types such as audio, video, subpicture aretime division multiplexed (TDM) within a common DVD program stream. Datafor each stream type is buffered before decoding by each of therespective stream type decoders. As illustrated in FIG. 8 d, thesebuffers can allow extreme variation in the time in which coded datacorresponding to one frame enters the buffer, and the time when it islater decoded and presented (e.g. to display). For purposes of buffermodeling, these stream types are deemed concurrent, although areactually serially multiplexed at the granularity of a stream packet. Ifa concurrent multiplex of the codebook would adversely affect otherconcurrent stream types (video, audio), thus leaving too little bits forother concurrent streams, the encoder may send the codebook far ahead intime during a less active period of the base layer.

Multiplex Method

The majority of DVD payload packets are consumed by a single MPEG-2System Program Stream comprising a multiplex of Packetized ElementaryStreams (PES) as depicted in FIG. 8 a. DVD packets (8004, 8006, 8008,8010, 8012, 8014, 8016, etc.) and are 2048 bytes long, but other non-DVDapplications to which PHD are applicable may have other fixed orvariable packet lengths. The flexible aspects of the of the DVD cell8002, 8102 structure (buffering, type order and frequency) aredetermined by the DVD author. The example cell 8002 demonstrates thedominance of video packets owing to the larger share of the bitstreamconsumed by video. The actual order of packet types within the stream isarbitrary, within the limitations of buffering prescribed by the DVDstandard and other standards incorporated by reference such as MPEG-2.Each concurrent data type within a DVD title is encapsulated in themultiplex as a separate PES. The program stream is an assembly ofinterleaved concurrent PES stream packets. The standard definition videosignal (packets 8006, 8008, 8016) is coded, as per DVD specification,with certain parameter restrictions on the MPEG-2 video tool andperformance combination well known as the “Main Profile@Main Level”(MP@ML). Other data types include Dolby AC-3 (8008), Sub-picture (8014),and navigation (8004) layers. Each PES stream is given unique identifierin the packet header. Room in the ID space was reserved for futurestream types to be uniquely identified through the RID (RegisteredStream ID) mechanism maintained by, for example, the SMPTE RegistrationAuthority (SMPTE-RA).

PHD would appear as an additional private stream type within themultiplex (FIG. 8 b), with an identifying RID. Because they appear as aprivate stream type, PHD packets can be ignored by older DVD playerswithout consequence to the reconstructed MP@ML base layer video. Othermultiplexing schemes such as MPEG-2 Transport Stream (TS), IETF RTP,TCP/IP, UDP, can be adapted to encapsulate PHD enhancement streamappropriate for each application. MPEG-2 TS, for example, are suited forbroadcast applications such as satellite, terrestrial, and digital cabletelevision, while RTP might be chosen for streaming over the internet ora Ethernet LAN. Program Streams are required by the DVD-Videospecification, whereas emerging DVD formats such as Blu-Ray have adoptedMPEG-2 Transport Streams as the multiplex format.

Codebooks are a significant portion of the PHD enhancement stream. A newcodebook or codebook update is optionally downloaded at the beginning ofeach scene. The other major portion of the enhancement stream compriseindices for coded enhancement blocks.

1. A circuit for enhancing picture quality of a video signal, saidcircuit comprising: a base decoder generating a base image of a standarddefinition picture; a temporal predictive interpolator coupled to saidbase decoder and generating an interpolated block; and a summing circuitcoupled to said temporal predictive interpolator, said summing circuitadding said interpolated block and a difference block.
 2. The circuit ofclaim 1 wherein said temporal predictive interpolator comprises acircuit for providing motion compensation.
 3. The circuit of claim 1wherein said interpolator is coupled to receive previously decodedpictures.
 4. The circuit of claim 3 wherein said previously decodedpictures comprise a first previously processed picture and a secondpreviously processed picture combined to generate a combined motioncompensated surface.
 5. The circuit of claim 4 further comprising asecond summing circuit coupled to receive an interpolated base frame andsaid combined motion compensated surface and to generate a temporallyinterpolated enhanced picture.
 6. A circuit for enhancing picturequality of a video signal, said circuit comprising: base decoder meansgenerating a base image of a standard definition picture; temporalpredictive interpolator means coupled to said base decoder means andgenerating an interpolated block; and summing circuit means coupled tosaid temporal predictive interpolator, said summing circuit adding saidinterpolated block and a difference block.